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(54) Abstract Title 

Correcting for distortion in seismic data due to ship motion 



(57) An apparatus and method for removing the 
distortion in marine seismic data resulting from the 
motion of the ship i.e. Doppler shift. The ship trails one or 
more seismic sources and receivers and moves forward at 
a known velocity. The seismic sources emit seismic waves 
that travel through the water and reflect off interfaces 
between rock formations below the ocean floor. The 
motion of the sources and receivers introduces distortion 
in the recorded seismic data that can be modeled using 
Doppler theory The data preferably is corrected for source 
motion independently from the correction for receiver 
motion. The seismic data is first corrected for receiver 
motion and then for source motion. The technique for 
correcting for source motion includes correlating the 
receiver-corrected data with a reference sweep signal, 
performing an F-K transform, performing an inverse F-K 
transform on a selected subset of the F-K transformed 
data, and computing appropriate correction filters for the 
data resulting from the inverse F-K transform. This process 
is repeated for all subsets of F-K transformed data and the 
resulting filtered data are summed together. 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 
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to a subsurface structure, the geometry or topography of the stracwe esa be <tamne& 
Certain topographical features are iodicstive of oil ami/or gas reservoirs. 

To determine the distance to subsnrfecQ structure «0, source 1 ?2 emits seismic waves 
115 which reflect off subsurface structure 60. The reflected waves are sensed (by receives* 
J 54. JSy determining the length of tarn© that the seismic waves 115 took to uzrsvel from source 
112 to subsurfee© structure 60 to resets 114, en estimate of the distance to saflbsudfsse 
stracture 60 can h© ofetsatnxsd. 

The receivers ussd an marine seismology are commonly inferred to as Cnydro^ooiss, or 



marine pressure pnomss, and are ui 

10 Synthetic piezoelectric materials, such as barium areoimte, barium aStooto, or 
outaaiobG^ am generality used. A sheet of piezoelectric material develops a voltage 
difference between opposite feces when subjected to mechanical bzMm&. Thin electroplating 
©a these surfaces allows m electrical connexion to feg 80 the d^ce so tto ite vollsage 
63A fee raessured. The voltage as proporttal to she ammmt offnecbsnkal bsmdieg w preiSHre 

15 chaste estpeirBeitsed by the receiver as resulting from seismic energy propagating through the 
water. Vmoas types of bydrophomss ore available such as disk hydropboates and eyBMrieal 



Two types of seismic sources are used ao generate seismic waves for the seismic 
measurements. The first sswce type comprises m impulsive source which genres a high 

20 easrgy, short tsrose duration) impulse The torn® teweeca emitting *te waps^ ®» 
and testing the reflected impulse hy a receiver is used to desesimine the distance to 
sstafise Btrocsare wsacfe investigation. The impulsive source end the associated 
Gspisiti©© BEd ps®ces&iss^ sysaem ear© relatively gampBe. However, the ESMsgaitude of eaeirgy 
irepared hy seiguMC tedmiqoei using tapulsave sources may, in some situations, fee hmsM to 

23 marine life an the immediate vicinity of source 1 12. 

The emdronmeatal concerns assocaased with impulsive sources has lead to th© uss of 
another type of seismic source which generates a Dower magnitude, wbmozy ®nsrgy. The 
measurement technique which us*s such a source as referred to as the high fidelity vibratory 
seismic CWVS") technique. Rather than imparting a high magnitude pressure pulse into the 

20 osean in s very short time period, vibratory soarees emit lower amplitude pressure waves over, 
a time period typically between 5 and 7 secoaads. Further, the frequency of the vibrating soaarce 
varies i&om 5 to 1 50 Hz, although the specifac low and high frequencies differ from system to 
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system Th$ frequency of the source assy vary libe&rly wife respect to time or noiB-liiaesdy. 
Tfeg &eqiu>es!icy variations are commonly called a "frequency sweep". The iStegiiectey swes^ b 
thus &s*w©gaa Sand 150 Hz sad S to? seconds in dwofiioa Tflk magnitude of Ifes eeisaiic wave 
osdUhatioos may vary or ranaiin at & @on$i&n8 sffraplofiude. Tfee smplitisde of the ©sdlatDonSp 
5 hwver B are much tower ihsn tfes magnitude of topuMve sows@s md $ms> te^ ©re fewer 
(^viiroinimanitoJ concerns wth the HFVS seismic tedinique. 

Seismie ships musl move fomard white seaside mmm&m®$l2> being ir©sord@d for 
msxsy o-esscsis. Referring sail! to Figure 1, the hydinsphones I H ©annexing mm> oM s&ress 
fiaeroubers provided on the streearcers ars plaesd inside a neoprene tebe (not showa in Hgwe 1) 
HO 4-3 isacto in diameter. The tube as then filled with sufficient Baghter-thsn-wsJcr liquid to tnaske 
Sfes s&resmsr Qgutzally buoyant. A le©d-in gecaioon 111 of Cbs streamer 1 10 apprcrasnMely 3C0 
feel tog o aasmbsr sireieh of aestiicOT approximately SO ©asters Boag trail Ib^Rfeera tfee 
eliip's s&eraa arcd the streamsr section 1 16 in which the receivers 1 14 are included. A deflector 
.pmvws 113 pulls the streamer section IJ6 out to sun approprialle operatiing wMtk Bteptlb 

33 eofflfroHersi (not shown) ore fastened to fdte streamer est various places ato^g its fleagtlk TTlfeese 
devi©3& g^roe ftfeg faydro&tsic pir<eg§uure ©axdl site bird wings 20 sfeat the Stow of water <ow^r tfesm 
raises or bwrs the stressa©- to the desired depsk ITha rilepsh that the Gonttofllesrs sedis to 
mai$&i® cos ba eoatroJlled by a signal sent through tEse sfcreOTor cabling and dans the d&jp&h eon 
lb e&iiKged ss desired. For the streamer's deptfo control §ymm to &roc8ion effeEtiwsBy,, ftike ship 

20 100 smsa travel fewas^ at. a sps§3 through the wstcgr of approximately four knots. 

Sseoffid, csreaflSKgs' 110 uaaolfly os a flwibOe sdbfe md tfimas die sMjp mniust smmre forward to 
fiaatasdlifli o dosistsd S?t@d 8epoffati©a towsa Cite so&rces streaan<ars, osd &@8w©@n fifes 
£3FesEfl£rg dJte^aivsB. "i^e spades : &ssto©3od so;u!u©ct aad sss^ameirs 5s 5mpOflt» m Efe© inaarfQ@ 
c^fiin25j>l<ogy a^d raaM mi vosy Wlaafe gdsmie mess«ss^(3a!l or© Baad& 

2S TMird, csoeinaflc shops oto deploy onultipte gtr@airors ming barvoaes that allow a f5u@d 

sepoinDtiTO to (be maitiiiained b^we^ra ^n^msrs. Ths&z ten/eus^ fore® tte si^eainissrs Bsa@rs!18y os 
ths kmt mxm$ forward. WitkmS efee b&rvsaes, ^h© s8r®am©re snay tesoma eateagled. The 
hm^ ^sed deteroliaes the amo^m of separatbn beiwe^ streamers. 

Foytfft, sdasaic ships mia^ cov©r as rmjcSi oosiiw swrfese as possible e&eh d©y becaam 

3© of sc^E of opsrstoBg the sMp. For these r^sons esw8 others, seasaiic ships must hksov^ 
fomorf wfealo teMng i^^surements @nd the forward sp©sd onust be reasonably constant. 
IVpltea! shop Bp©sd is approximately 2-3 meters per second. Because the streamer is deployed 
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behind the ship, the acaree and receivers si so move at approximately 2.S metes per second. 

AWfcougb ship motion io necessary as described above, the motion distostD or ''smears" 
the ©squired seismic data. Broadly, smearing results from the feet that the slluip, ami thus the 
scum® and receivers, move whale data collection takes place, lit is generally recognised that 

3 the smearing etffect of ship motion en seismic data results from two analytically separate 
phenomena— source motion and receiver motion. Although the receivers and source are pulled 
behind she ship and thus move sit the same speed as the ship, she effect of scares ostction on the 
tea is usually aaalyzed independently from she effect of receiver motion. Source motion is 
less ©fa concern than receiver rotation in impulsive source-based seismic systems because the 

10 seance EEOves a negligible amount during the brief impulse emitted fey the soaree. Data 
smearing in a HFVS system includes significant contributions from both receive? and scarce 
motcon. Thus, the HFVS-acquired data Ami id be corrected for both (receiver and source 

SfcOBOOaL 

TBte high costs associated with operating a seismic ship require that the ffliethods sad 
J5 prosssures used be efficient It os thus desirable to maaimise data collection ia as short a time 
as possible. Bosaase of the length of the frequency sweep (typically 5 seconds or more), 
HFVS sources are typically activated ©very 10 to 20 seeoisds. Because of the ship's speed (2-3 
meters per geeoasd), a HFVS source onuss be activated mo sooner than every 25 to 37.5 rosters. 
Although more data fai oose location could be acquired if the ship were So travel as s slower 
20 spssd, streamer control would be lost and less ocean surface would be covered each day, 
thereby jncfeasimg she cost required to make seismic measurements of a desired section of tlhe 
©esoa. 

sssosdba data. M an aracte eaatitled "The Efifccls of Souse© and deceive? Motto ®a Seagate 
23 Data,™ by Hampsora and Jafcubowicz, Geophysical Prospecting, 1995, p. 221=244, a method for 
canTesJj^g few receiver and source motion is disclosed. Although the method of Bampson sand 
JataabowicE saay have theoretical merit, the method is impractical for use with conventional 
marine ceismic systems as it requires the H1FVS source to be activated with a temporal sntd 
spatial spacing that is impractical. It is well known that for a wave traveling with a velocity V 
30 tGarough a medium such as water and with a frequency of F (Le., the number of complete cycles 
of the waveform per second), the velocity V is related to the frequency F by the length of the 
wave, relfened to as the wavelength (X). The relationship is: 
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V^IFoL (1) 
Titos, the wsvdeHngfth & as V/3F. In wer seisnaic waves propagate with a Etrcowri vdotity of 
^proramsfcely 1500 sneers per second (approximately 3325 miles per hour). If tbe higher 
fr©qaeaey an a cta/sep 5s assumed to bs 60 cycles per saeoad (or 60 "Hz"), fl&e ^veleragth of 
sueb a seismic is 25 m&m (85©0/6O). To avoid & certain typs of date distortion to^wa 
as "oliasiiag™ tfl&s source ram be activated 08 s spsdag of at leas8 one half of ttfes ^aveBe^lh. 
Tkm^ for Hampson acid Jakubowicz's method to work tfee viferatoiy smirees smusa fee stitfosted 
at <srvery 12.5 racers, and preferably sooner. To activate a source 81 §aeh 
spodags, ship must travel much slower fth&m ots preferred 3«4 rovers per second!. 

it wodd be advantageous to provide a practical seismic system for use in sn&rine 
sppBieotioias to ean corred sbe dsta for' the motion of the ship wSalbmst tte defieteiaeiies 
ioksnm in the Jfeuropson am! Jakabowfe method. Such a sysienfo preferably would ©orsrset for 
b®th jrecdves* and source motions aacfl do in a eost effective maniaer. IDespite the ^parant 
edTOi^teg^ to data all attempts of cievelopiftg suefo a sys&<£m teve jailed. 

Tlfeg pzobfleims outlined ®bw© are in terge psrt solved by the seismic ratsasuresnenK sotd 
processing system of the preset invention. Tte seismic raiig&SOTemsEit and processing system " 
disetessdl bea^to removes the distortion in marine seismic date resulting fiosn ftfee motion of ttes 
cMp. Assorting to Sfoe invention, the ship Qows W&and it one or more seiaraie £OTT&g£ Qiad 
stre@giB@ffS it saoves fomsni at s constant veBodty. The sstesaie sources mk seismic wws 
ftei travel flkough ftfoe water su&d rfea off i?nterfkces be(twe©n rock formations bekw She 
SHcot. lie motion of tine sarogs &®g revivers imrodtaice distortion iso fbe r@6ord@di 
estate dksa tfest ©ar& be mod@ted nasing Bopplsr th@©sy. Tte date preferably is scraped! fis? 
erar©G i^ffifiiom w^pooAos^ &®m ofe© ©onre^ooa f®r r®s&v& snoEoon. AssjmSifig uo Cte@ 
psTsfeiisd embodiment, tfes skgisraie dsoa is conm^d firss for reedv^ motion wmg any of a 
vmety of tedBn5qt&^ oad ttesn for soaarce fiacSDW. 

Tlte t@cWdp© for wwcdoBg for sewse mnotiow includes correkiSng the r©P©5v©r- 
co^o^d dfejto with o Deference sweep signal,, performing a fcrsnsforcn (mdi m m F-K 
tmrforai). p^fonnnisig an inverse transform ($meh ss an inverse F-BC transform) odd a selected 
edb&§8 of tho traimsformsd data, and computing appropriate corredioon filters for the dlsta 
r^teSpg i!h)fin nbe iiirras© F-K transform. The inverse transformed data Gonresponds to seismic 
easrgy tfeaa travails upward from subsurface structures at a particular angle r&fenred to s dip 
oagSe. Appropriate DoppSer eorrteoion filters ore computed for each set of inverse transform^ 
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data and the process is repeated for all subsets of F-K transformed data The Doppler filters 
are applied to the seismic data, and the filtered data are summed together. 

These and other advantages of the present invention wfli be apparent to one skilled in 
the ait upon reading the following detail description of the invention. 
5 A better understanding of the present invention can be obtained when the following 

detailed description of the preferred embodiment is considered in conjunction with the 
following drawings, in which: 

Figure 1 shows a ship for making seismic measurements with a towed streamer array 
including a seismic source and multiple receivers; 
10 Figure 2 shows a seismic measurement system in accordance with the preferred 

embodiment of the present invention; 

Figure 3 shows the preferred method of correcting seismic data for the distortion 
caused by motion of the source and receivers; 

Figure 4 shows exemplary pressure data from multiple receivers and the distorting 

15 effect of receiver motion on the data; 

Figure 5 shows a preferred method for correcting seismic data for the distortion caused 

by receiver motion; 

Figure 6 shows an exemplary plot of seismic data in the F-K domain; 

Figure 7 shows exemplary shot records of multiple receivers in which only data at a 
20 constant dip angle is included in the shot records; 

Figure 8 shows the relationship between apparent wave velocity and true wave 
velocity; 

Figure 9 shows the preferred method of constructing and applying Doppler shift filters 

to the shot records of Figure 6; 
25 Figure 10 shows the geometry associated with a moving source, a single point 

diffiactor, and a receiver for calculating the amount of Doppler shift caused by the moving 
source; and 

Figure 1 1 shows the preferred method for computing Doppler shift corrective filters for 
multiple diffract ors. 

30 The present invention corrects seismic data collected by a marine seismic system ft>r 

the motion of the towed seismic receivers and sources. For simplicity, the technique will be 
described with reference to a "diffractor" (also called a "scatterer**) which is a reflecting point 
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located at the physical interface beautfaen contiguous subsufffoce formtioas. Because a 
subsurface interface m comprised off many point dffirsctor^ the entire interface <5©n be mapped 
by finely superimposing tfeg results from eaoh poiat daSrcsaor. 

Rgfemag j&ow to Figus© 2, a seismk system §0 comsanaGed in accordance woth the 
preforad esribodinaent generally indues & seiismie roe&sureffifteot arad prosesssiag system 51 1 q 
user ais^ua device 59 (pgferobJy a keyboard, button switches, aiad cosufifd hacks), a displtoy 
dovi©3 32. ©R3 or caor© seiss&ic sowcss 1 12, ©nd outs or more cables (aOs© colled mmmr^) of 
s&imte receivers 11 14. A streamer of receiver is also referred to as a "gather" Tfcs seissate 
measure**^ ostd processing system SI includes a prcee^irag mk 53 coupled so a dusa storage 
unit $4, o&td csuaras w& recdver Interfuse mk 36* It shcwdd be recognised thst fcSne setaife 
meaDuresftgai oyssem 30 may imcilttde other cora^Qeffits not showsu is Figure 2. Tfe© user iapfl 
devke 59 allows a user to input commands and conSguimtijwm ?nforiraia&5offl inao the cy&em SO. 
The dicpHay devise 32 provide wual jngpiresgm&ti m of daftn, Gonfigurstion itsfoirEBrfofl, oad 
status informatics to Eke uca. The source md receiver preferably eoyple to ahe seismic 
measurement system vis Sher-optic cables S7. The sowe 112 includes my mtoM® setsffiate 
scunee such as HFVS soosrees sad impulsive .wrass. IRese&vers 114 telicde§ suitable 
hydrophone receivers including piezoelecforic-based devioes Off any other gyiiiteble type of 
sdssaie reiver. 

The processing wait preferably controls the option of the seisinak Efi^w<»eM 
system 50, stariiag dstn in storage wild 54 (which pireferably as a ragrtetic (tape, § ted dis?L ®r 
CD ROM drive), and cofmroSlorog flhe opsmion of scots© 112 sired receiver 114. Sdssafo 
signals detested by the receivers m trsasmfcted to nine sesgmie onDea&jremesKt sysSem, jprosessdl 
by ^no^booag 52 asbd'eittrsd ia Gtorsg© eaait 54. 

S^ferimig sow to Figures 2 oad 3, artd 05jploiRj©d! in mraore detail Sa 4fe@ di2SQ3Ssifi)Sii tefl 
Ifeltewa, tks gsasmic meanw^eflil and processing sy-ssem preferably g®tt®&& tfc® record©^) 
sd^io dMa for ghe snotiom of the receivers 114 and idfcs smsree 112 scoordi^ tthe 
isa©llhodology MJtastrsted om iow ckirt ISO. Altem2£svely B the seismic data esua be sSered @s& 
ssssgBixgtBc tape or dick aiad tonsfeired to anc^bsr computer system (Tor aroalysas accordifiag to fife© 
fceestags of nh© pxnsfennsd embodiment ot & loc^bn remote i&om the sdMaic ship. Tibs 
prefer dats oorr©5t5ora niethod corrects first for the effect of receiver motion in step 160, mA 
Uto art data to the effect of source mofiion in siqps 170-240. Each of aJtese s^eps is 
©spBoxaed hdom. 



F-7970 



8 

^<a)DT^eft5^ca far I Recciiveir MoCioin) (ste p 1160) 

Referring now to Figure 4, a gsafeer 1 14 of receivers 125, 126, 127, 128 ds shown with a 
pressure signal 120 recorded by each receiver. The pressure signals 120 are referred to 
coHeetovely as a "shot record/* Time is represented along the vertical axis and distance is 
3 represented across the horizontal axis. An exemplary trace is shown for one receiver 125 ©nd, 
for simplicity, a straight Ibe is used to represent the remaining shot records. 

Hf the receivers dad not move while recording the shot records, the traces 120 would! be 
recorded at a fmed location and therefore would! be a junction ondy of time, sad stof space. 
Because the receivers are towed behind a moving ship (assumed to be snovarag to the right in 

10 Figure 4), each shot record is recorded as a fonctoon, not only with respect to time, feflt also 
space, as indicated by ftracss 122 for each receiver. Traces 122 represent fcraoss 120 as the 
reeeiwar is pulled behind ahe ship. Thus, each data point om the shot records 122 repress Ohe 
seismic pressure signals sensed by the receiver at a particular point in time and §pa©s. 

IRefemng still to Figure 4, each receiver is assumed to be located tt [posfciota cq when 

15 the shot record begoras. Thus, receiver 125 begins at location) H2% Receiver 126 bsgins at 
les-jjtooa rllKs. receiver 127 at locatson r!2?0o and receiver 128 at locatooa jtU2% The 
distance between the initial Bocation m and the ending position is a function of the speed of the 
recdivere. It is assumsdl for purposes of this discussion that the speed of the ireagiveys, as well 
as Uhe speed of tifee sources* is the saime as the spe©d of the shop, although in tlfceosy sligbt 

20 differences Sa the speeds m®y exist due to such factors as the elasticity of the streamer 1 10. 

Sfoofc irecords 122 ore repmssarted in Fipure 4 as straight diagonal Icoes. The fees 
(repr©E3ntaag pressure waveforms) are straight b^causs the receiver spesdl is assumed! to be 
Kwtoat. If (fee needver speed 5s op, then the position of each receiver at any tiro 8 <&ms§ a 
shot resold 5s m + od*- T5he linesr slanting of the shot {records 122 is equivolent to & time* 

25 woont ^psUal shift. If u<%«qQ«S(t)jrt%»8) represents the magnitude (pressure) p> of the shot 
record os a fijaciion of source Jobation a, source speed seismic signal Stft) produced by ahe 
source, receiver location or, receiver speed nr« snd time 8, then time variant spatial shift csn be 
HBBiiaheraaticffllily modelled as the convolution of oK&9ros^t),r f niftlt) with a "Dirac" deita tfuaefcion 
(also referred to as a "unit impulse 9 *): 

where the ° operator denotes convolution and 6 denotes a delta function. The convolution of 
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tem fcrsefcions (s faction presents a series of values at various points in time or space) is a 
fatam maGheonatDcal operation which involves replacing each element of one ftanciion with on 
owtpaut fiutetfon seofed according to the magnitude of she input ©lenient, and then 
superimposing the output vote. For a more detailed explanation of convocation, referanee 
(saa fee made to "Exploration Seismology/ 1 by Sheriff and Geldurt, published fey the Press 
§y©di©aft© oftho University of Camlbridge, 1995, p. 279-81. 

nice spatial shift represented by 8(tffl> + oM) to equation (2) can be ^amoved by 
convolving the result in equation (2) with a spatial shift in the opposite dfasetoon. The 
cofyectioro for receiver raiodon h therefor©: 



la eqiesirioia (3), eanvoluftbm of the spatially shifted shot record with the delta toetson 
vtft) iresuto in o shot raeond kid Ghe ir<&cdiver teen stationary (cv=0) all position sty. thus, t&e 
effect of receiver motion on the shot record 5$ neutralized hy convolving the shot record wiofc a 
deka fbn^fon representing a spatial shift. It should fes recognized thai the foregoing aiaallysis 
involves tfwetioas mathematical operations that occar as tactions of ttixne wd spaee (te 
so eall&IL ttscsm© and space domains). 

Otto my§ to correct the shot records for swiva motiosu sire available. For OTsannple, 
She conmion provided! on equation (3) earn alio be represented in the §eq«cy domain in 
vM$h all dunctiotaa vary wfch ffequeifocy, not time. Functions can be converted torn theiir tonne 
osd sp&©3 dosjBasn representations So flhe frequency domain using a snatheraafel operation 
.©fllsd a transferal. Tte (fhaqpgnctes involved with aaeh Fourier graasfoms fehDde 

tcaaporaB and spatial frequences. The Fourier transform of the deha function, 8{ir® ° fe <2 

'® B * D » fl ^fae 8 repes^Jts she onaaginaffy numb®' (the square root of -IX & represses She 
spatial iErsquency (also irafemred to as oh§ wavraamtotr) and ^ 5s a known constant. Jt is well 
known fhtt convoJallKMii in the time and space domains is equivalent to multiplication in flhe 
frequency down. Tims, the spatial shift introduced in equation (3) to counterbalance tihe 
spatial shaft eaiG©d by the receiver motion can be represented in the frequency donasin as ftbe 

product of the Fourier transforms ofthe shot record and 
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Where fftf, h) is the Fourier uraaisfoirra of the sbd record and is a function of tempoiral 
JrhsqpueKgy toad spstSal frequency 2s. The symbol "o* denote multiplication. 

Rafefriag Km to Hgaares 2, 3 end 5* the seismic sness&remeat system 50 sesaovea the 
©Seen of receavar fflootiosi uisSsig equation (4) by first computing the Fcwr transft&ra of the 
s sfoot records in sfcep 162. The sdssuc measurement and processing $ymm SI eoffiigutes the 
Fourier transferal) laslng aay « of a variety ©ff knows techniques suefo as the Fo& Fourier 
Tiramfom la tftep 164, the seismic ffi&ei^ei?^ aiajd processing sy$t@fift SI mosihtipH^ the 
Fourier transform of the shot records by the Fowier transferal of the delta tosfcitaa of ttpttioa 

(3) represented es e^ KDWJ . Finally, in step 166, the product from step 164 h ooavesrssS 
1* badk 'mo the tirr&e and space domsh through m operation) referred to es the tore Fc&rter 
iimmfom which ol©o is a Ecnowm t©ehaiqu@. 

Another amethod for ocsrecaing for receiver rotation 5s described with release to ¥igm& 
4 to eoarect She shot mssord for the inaction of the receivers. This roe&hod will be descoiibed with 
irefereneg to osae such receiver, such as receiver 127, In this method, the ssksftie 

13 processing sys&sm SI selee&s dista from a reosiver white the receiver is in^sr tics location at 
^dh ahe shot *e@ord 5s to be feed. IT© fbrth@ to record) for iocaiio© crlL27<fc for ©fcasa*|p1e 0 the 
ssis&tDie measurement system selects the portion of shot records &obssi raseavers 127 0 H26, and 
125 when each ireseiver is near location ri27(|). The portion of the shot records to lbs selected 
by seismic messyresaeaa andl prao@ssii?g system SI is idesjrffned by reference sarafes 82?©, 

20 B27b 0 aoBdl 127c. Tfeas, seiaimfc meMimmeiit system selects the initial portion 127a of she flfeot 
jrasord tos iraseforar 127 uotil that receiver moves a distance ^proximately ©pal to <M@-Wf 
afcs mmy location trE27@. At than poooat, sefeonte ffifi&granM^ 

pios^ftlg systesaS 11 geSeetethessiaddlepoftflOfii 327hof8he^trecorfM^fit5S^w^ 126us^Sl 
to reaver oko eeomss <Mse-Hf th® gnrop imepval away from location rS27§. FtnaSBy, the 

25 lost portbra !27c of the shot record from receiver 125 is selected by the sys&ssa 51 . 

ITfee H©gtfo^te described ak>ve ®re caomplery only of the flDiesbods for ©oirrecttfling for 
recover motim wd the iaventooa is not mt ended to be Jimited to my prntimlm' method. 
Prdk&% folkmhfag recdver motm correction, the seismic measurement antd processat^g 
sysfcean 5S eorrecSs lifee dfga for the sotsrce motion. 

30 QpFTCCftfaifl for S^natff5Q M ojimigteias »7<^24<ft 

Refeoiag to Fij^ar© 3„ in she preferred method 150 for correction of the data for soiarce 
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motion, the* seismic measurement and processing system 51 comsJattes the dm (now corrected 
for receive snotioa per step 160)w«hAeHF^S referent swe^ sigpaJ. Th© HFVS r©fe^© 
awssp sigwl casa be aay sweep signal desired ausd snay include linear frequency sweeps 
(frequency changes at a ©ooistant Me dunsBg the sweep) or non-liiaear froqugtigy sweeps 
5 (SrequOTGy changes ®t a varying rate during £hs sweep). As erapllaiaed below, the eorretetioisi 
step HO is aeeess&ry m & HFVS system to compress the retetiively Jong swesp to a short 
duration <3voiai. 

Tfee earth csn be thought of as a filter of seismic energy. That is, if seismic energy is 
input into the ©arafo, a receiver portioned ©a ihe sarfeee of the girth will receive setoic 
BO enessy whose character has been altered by the earth. The various Gaetors that ms&ify tfes 
seismic wove as 54 pasess through Hilts eauth include: 

(a) ihe iae&r ihe seuree where the Presses md absorption of <m<&gy ohm 

(b) the iresponse of the dimeters comprising the subsurface interfaces (the signs! ttefc 
BS ssSsmie worfc is sostemted to fiadl); 

(c) the neor-surfsce sorae, which has a disproportions effect in modifying the twiv@; 

QSDd 

(<dl) ©dJdfekjrf sacdifymg etffem because of absorption, w coftvsrsDOia, smkoplos, 
ascd di$&©ctio@s, m& nhe BEse. 
20 Is pssafee, the ireeeiwrs ireeord ksj only prfmairy seismic reflee&iiofls, hut also multiples, 
«Sa®®§toom, scattered wives, etsflect@d refractions, surfece waves, and she Jitoa, all werilappifiig 
itatfcige. 

Ganssdily, a {Elfcer is a sy&em that prodbces an cudpua sijgnal'fixf £ given inajpaBt sgsaal. 
ISds @&Jj?y3 sigpal G§a hs caHealsK^d if 'the IkjmOse response for the filter fe fomm. Tib 

2S tapifce msspoo&Ds is the onSpMt signal proceed by the filter for © given taapulse iapaat sigrool. 
The ©aatpast sognai as simply the input signal convolved with the impulse rcsjpon&ge of 6h© fik&r. 

TDfee seismic signal detected by the (receivers represents fihe input reference signal 
iftffcesacsd by the iteors described above. Seisimic data (or "seismogr&nas") is ©sefful to 
<&&gmme the location of oil and gas reseavoirs when the dm represents the 5a^a8 ref^omce 

20 sogE&B ccted Bspoa only by the difftaetoirs coimprisiflg the subsurface interfaces, as contmgied 
with a ref@ngK©e input siga&J thae k also innuenced by the above-described sigf&al altering 
SjcSkts. The etffecft that the diffraaors hsve on the stisrme waves propagating thirou^i the earth 
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is refers to as She impose response of she earth. Because of the Atonal signal altering 
faeaors described above, the seismic signal received by the receivers in a HFVS bsara Me 
reggsBsbSsaos to the impulse responds of the earth. Seismic work is intended to determine the 
impuBse response of the earth, and thereby remove She influences on the data to are act of 

5 icsterest to seissnologte. 

To fesjiove the long sweep station front the recorded data, the seasme m^mmma 
system SO preferably sorreSates the resordsd data with the reference sweep signal. The 
©ontslatiioa of two data sets 5s a known mathematical operation to which on© data cet is 
displace by varying amounts relative so the other data set and corresponding valtags of she tow© 
10 seas are multiplied together and She products summed So give the valve of She correlation. In 
Step 170. shown in Figure 3, the data from ssep 160 that has been corrected for receiver motion 
is correlated wish the refesetac® sweep signal. 

In step ISO, an 1F-EC transform (F refers to temporal frequency and EC refers So spatial 
fequeoDcy or waveaumfcer) as performed on the cmrdated data torn step 170, altagh other 
13 suitable toansforms, such as the Laplace transform, radon transform and %-p transform, can 
alto be ucsd. The F-& transform is a doahBe Fourier transform in which a signal that is a 
Enaction of time; % and space, a. is transformed so a signal that is a Sanction of frequency, ft 
and renumber, Es. lbs (transformed signal! can be plotted on a plot referred to as an F-K plot, 
gyieh as afoat shown in Figure 6. Converting a ftincrion from the Sime and space domain inso tfije 
20 fcequ©®ey wti wovesMa&er teaaifl k referred! so as a forward F-K transform. By analogy, 
©saraersing o tetion from the iequency and wavemamber domain hack too Sfe® Sim© and 
^oe§ domain as reltad so as an inverse FX transform, the forward F-K ussrform is 

i a double integral! as: 



23 wfliere IFlJCssfi) is she F-K srarosforra of The inverse F=K transform (pefffforsned in «ep 
300) is represented as: 

Referring again to Figure 3. in step 190. the seismic measurement system seleols a 
e@3i8tarot stoae dip slice of daaa (described below) from the F-K plot. This ssep is best 
30 understood with reference to figures 6, 7, and 8. Figure 6 shows an F-K plot of a sransforatsd 
afeos record font Figure 7. Frequency measured in cycles per second or "Hertz" (Hz) os 
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presented oa (fits vwd axis and wgvemunber measured in cycles pgr i^er id irepireseflted 
oa life ton^mall oris. The F-EC tc^otrmsd doaa os represented fey porate 191 5a the F-& 
plot 

Evesy staiglht Bne, saeh as liaes 11 94, 19S, 196, begigmiimg Ifrofifl the erigja of the F-K 
3 axes wd exitendiqg eartwtrd represents seismic data vMh a pastJcuilsir apparaft ve!odfty. 
Fujrtfesr* the slops of ©as&a such stoigk lixae 5s ©spaS to on appareatt veltadty. Refkraag to 
F5$gre 7 0 insseSvere 325, 1126, 127, J 28 are skovm wM* a seismic w*re J 32 |po^pga&ig through 
the <§orth (Mudi@g waiter) in (fee direction of asrow 129. Litfte 130 represses the direction of 
propagation of seismic wave 132 aad forms m wish vertical line 834. Ttot ssgl© is 
10 referred to as she angle of approach, apparent dip angle, or samply dip angle sad is denoted in 
IFigptgs 7 estd 8 as ©MP. Doe 130 thus 5s referred-to as the dop line or Dime of appffoach fer 

piarposes of this sppffieatioo. 

St@f©riqg to Fugwz 8, straight tone 133 is perpendicular to dip ifo© 130 asssi repressors 
ceheraasiieally th& ww&oRt ofwa 132 as fcy smvel upward at oh© dip aagi® ®MHP. 'Site 
23 wi^sfiroa? 133 propagaaes sp through she earth with a certain veJociity rdfarod to as she tone 
vetoc&y, V^e. Tfee toe velocity of seismic waves propagating through wmer as 
approximately 1500 meiers pgr secomi (3325 miles psr haw\ in general is ooasiden^l to 
be o WW. TFsro wkxfom e&a be easily d&eroined usbg g^y orq off a vsrfeiy of taaowsi 

3fl> Bsfemag to Figwe 8, the kmzontal cosnipoinm of the tnae vetedty vector fa referred to 

o§ cfeg opparea* vdoday, Vqpp. The apparent veibeisy, fc: 

f ■:. .-- ■ Vo^ip ^ V®me /pm$mspl>° ? .. CO 

Wl^^ia w ii§tfce^ Tfoe apparel vdo^haspfa^<^4^fi^^ 

to 8teS ft is tie velocity of the seismic wave 832 as detected (by she reedms. As mvefrant 
S33 Kcoves upward, receiver 128 will delect the wsveifrom before recdver 227 deseed ft. 
Fmsr^^, hscesase of the distance (beawesn receivers 827 imd S28 and tonwe iwwval heewaen 
wl^n (he wavefr(« os dexeoed hy receiver 128 then receiver 127, the wovofroimt will appear to 
hs towKiag hosisoataSly with velocity V 0 p>p. 

As &m be ssen by ^uadion (7), V^pp is onverse!y pfopon5oaal to She sine of the dip 
so angHe 0©iup, givea that is a constant. Thass, each straight Bine in the F-K pbt of Figure 
6, elte slo^e of ^Mch 5s defines & dip angle, ©©up in Figure 7 and 8. Moreover, data in 
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the F»& pfiot of Figure 6 along a straight line, such as liae 195, represeats only the seismic 
eaagy that propagated up through the earth at a particular dip angle, &d excludes seismic 
gagfgy propagating upwards at all other dip angles. 

Referring now to Figures 3, 6, and 7, the seismic measurement and pracssaimg system 
5 5 1 preferably corrects the data for source ratoiioa by selecting a constant iisae dip slice of data 
torn the F»K dotrcain in §aep 190 (Figure 3). An exemplary constant tirm dip slice is atowa fa 
Figure 6 as the portion of date 191 bounded by straight lines 194 and 196. Because to 194 
aatd 196 define a pae-ahapsd wedge in the F=K plot, the data comainsd tetween liaaes 194, 196 
is [referred to as a constant tisn© dip slice or pie slice. By selecting a pie slice of F»K tea and 

10 inverse F-K transforming she selected pie slice data an step 200. the seismic messa^^t and 
processing system 511 selects only the seismic energy that propagates upward through the earth 
within a range of dap angles deemed by the slopes of lines 194 and 1% Thus, according to the 
preferred embodiment of the invention, a constant rimta dip slice of F-K data ss selected in step 
m astd inverse F-K transformed in step 200. TCte size of the pi® slice can fee est to whatever 

S3 saas is desired and as generally a Junction of the accuracy desired. The saz® of the poe slice thus 
relates to a trang® of dip angles, ®®m? +M>M\?. 

Tfee resaah of step 200 as a shot record that has been corrected for receiver oraotaon and 
that represents the seismic energy that corresponding to a range of dip angles ±fi&WP 
that are related as described above to the- apparent velocity defamed by the pie slice. It should 

20 &s c^sagaksd tfegt the smSc eaeany at dip angle ©MP includes a nperpoatioB of seiSBSaie 
wives km iretasd off millooas of difSraetors along Sine B0. Usimg pmeipfles growded 
m classical! Bcppte absosry,, the data can lbs corrected for source motion. 

T© BSBdsrgSad th© wpgBMktti Boppler theory, referents© is ma^ to Fapasti 7 m 
which a sow© 1 12 en@ves fiom [location ssatt&e beginning of the HFVS fiequefieyswaep to 

25 location st the end of the frequency sweep. Point diffraetors 140, 142, 144 represeaat 
esemptey difSractor locations along line 130. Lines 145 arid 146 repress^ the direction 
ssisaate waves travel frosn the initial source location and the ending source .location Sej, 
respectively, to point diffractor 140. Similar lines can be drawn for seasanie waves traveling to 
dufractoss 142, 144. The seismic waves reBecied by diffractors 140, 142, 144 travel aapwaid 

30 along Dine 1 30 with dip angle ©imp. 
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As shown, the source 1 12 moves from left to right and ftbaxs moves away from diffroetor 
140. Because the scarce moves sway from the doSroetoir, the ipericd of the emitted iSrequeotsy 
mssp source signal will appear longer. Akematively, the length of the Sr©<p@n>ey w§f will 
spp@s\r Go hs longer firor&D the vantage point of dfflfi&etor 140. This change to flrequefisy and 
length of the frequency sweep 5s referred to as frequency shift winder Dopier theasry. to tMs 
©sai/aple, feowever, the scares approaches dififiraeier 546 Glaring fifes fie^uensy sweep, £®d to 
t&3 fespeflcy sweep be©oo*ie$ shorter from the vantage poiai of dif&aeto? 344. SMsrastflff 142 
is befow the midpoint of the source's trajecimy as St moves during the froqiseffDcy sweep, sad 
tos flfeer® it zero ®& frequency shift associated wi&h ditfFracftor H42. Moreover, the sSstostioso 
dbe to soasirce motion cara be represented by the magnitude of the frequency shift isstog EtoppSer 
theory. As will (be seen bebw, the magnitude of the Boppier shaft caa be computed for eodb 
difiteesor Doeatcon, or raige of diffmctor locaaioas, &rtd appropriate filters ean fee designed to 
softest tihe daia fcrthe dSstotftiion. 

Referring ms>w flo Fogaare 9, the preferred steps 210 to Gostmpute m& apply Bc^pJe? 
©0fir@stl©ro filters to ©orre& for soisrce rattSoin Sractodes first computooig the Mgotede of the 
disfcetftSom to step 2X2. Because the Doppier frequaeflcy shift alters 8be fengJh of the frequency 
sweep g& ©sdfo ditf&Bceor location, the ffnagnitede of the dnstortBOBD dose to sows© motion eon 4 be 
g^r©g@nited fey ©o^atixng the change in the length of the frequency sweep for each diffltector. 
He damage in k^gtih in the fr&juency sweep, unearned! m waits of milliseconds, 5s instead to 
os dfflsftOT (or coi^essifofl) §imd thssE, in step 212 the dilation 5s computed (for eaeh diffEBrast©^ 
The <$totfea wsfes w&h diffcetosr location aid is timss divided oroto time gsaes bub step 214 es> 
ifaft tibe E&fcfiate measurement osad processing system 5 1 can provide o conr@g&» filter for @oe& 
flosss j^fie fa step 216. FteiSIy, in $ejp 218 ths seismic measrosinient qom! pmssgsSiag systesaa SI 
z^^fiigs tibs ecMvecftoGsi SSters So tte afeot reoomd to eon®es fi«r d3ias5ofii. 

E^toang jftaw to Figure BO, the geometry associated with dfesrivbg tfes dildtiw fesr a 
diSlresSor J40 inxskdes soiree 1 12 caoviag &om ®n imfal lotion at fthe b^gimirog of a 
Snetpeaey s^e^) to eadkg position at the end of a frequency sw©gp. BMf&actor 140 is 
lo^t@d st a d<^&) Z (be!ow the steiooaary ireoeiveif 527 md distance % sway from the ireseiver. 
The &imM& 1 irepinesgms ^he distance between the receiver sad the source 112 at its initial 
jposWoft c®. The sagle as the angle from dip line 130 to the horizontal mis. Angle 8q» as 
ff^at^ to the dip angle ©©jp as 8p ° 2$ - $hpdp. Thus, once a constant pi(g siice of data is 
geleeied from ^he F«K domain, 0 r is defined. Seismic waves from the scarce at location s$ 
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travel along line 145 in the direction indicated, reflect off diffractor 140 and travel along line 
HO to a receiver 127. Similarly, a seismic wave emitted by the source 1 12 at position t~d 
travel along line 146 in the direction indicated, and reflect off diffractor 140 and also travel 
along line 130 to receiver 127. The amount of dilation is calculated as the difference between 

5 the time a seismic wave takes to travel from the source 112 at ha initial position s* to receiver 
127 and the time a wave takes to travel from the source to the receiver when the source is at its 
ending position w Recognizing that the seismic waves take the same amount of time to 
travel along line 130 between the diffractor 140 and receiver 127, the dilation is simply the 
difference in time a wave takes to travel from the source 1 12 at location so along line 145 to the 

10 diffractor 140, and the time of travel from the source at location a«j to the diffractor along line 
146. If T* represents the former the time along line 145 and represents the time along 
line 146, then the dilation is: 

D1L * Tmj - 1* <W 

where ML is the amount of dilation. The dilation value DEL thus is positive when TU* is 
15 greater than T* when the source is moving away from the diffractor) and negative when 
Tm»j is less than 1* (source is moving toward the diffractor). 

Referring still to Figure 1 0, and applying the Pythagorean theorem: 

20 and Z^XuatO,) 

where V is the propagation velocity of seismic waves in water (1 500 meters/second), T is the 
time a seismic wave takes to travel from the source at sO along line 145 to diffractor 140 and 
along dip line 130 to receiver 127. Equation (9) can be rewritten as a quadratic equation and 
thus can be solved for X: 
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where: 

-B±jB i -4AC (U) 
ZA 



30 
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C-H , +V 4 f-2H*V , f ( i2) 



5 

B'^HfV'f-H 1 ) (14) 

It should be noted that X in equation (11) is a positive value if the dip is down and to the light 
of a negative value if the dip is up and to the right. Using equation ( 11), X can be calculated at 
time T, and once X is known. Z can be calculated using equation (10). Calculating X and Z 
10 for each constant 6 r provides the location of a diftractor tor a seismic event at tune T. 
Equation (8) can be rewritten as: 

(15) 

DIL * J(H +«,J a - X) 3 +2 l - J{H +Z ' 
where u. is the source velocity and T«. is the time length of the frequency sweep. As can be 
15 seen by examination of equations (1 1H 15 X the dilation DIL is a function of the location of a 
diftractor (X and Z). T, boat speed u» and the length of the frequency sweep Th, Further, 
equation (14) only has a solution when T > H/V. 

Referring now to Figures 9 and 11, a plot of dilation DIL as a function of time is shown 
for a source 1 12 and a receiver 127. As shown, positive DIL values are plotted to the right of 

# 

20 axis 188 with negative DEL values to the left of the axis. Although, correction filters can be 
calculated on a sample-by-saraple baas (which may be preferred), satisfactory results can be . 
achieved in leas time if the dilation curve is divided into segments and correction filters are 
c onstru c ted for each segment, rather than for the samples that comprise the segm en ts . Thus, 
according to the preferred embodiment. dilation DEL is discretued into segments 252. 253, 

23 254, 255 and associated time gates 256, 257, 258, 259, 260. The size of segments 252-255, 
and thus the size of time gates 256-260, can be set to any desired size. 

Numerous techniques are available to construct appropriate filters to compensate the 
shot records for the amount of dilation in each time gate. For example, the sweep can be 
resampled to At' where 

30 T 

M= — a — A/ 

(16) 
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where AC as iBne sampling period for the shot record. After resampling, ih® » sample rase 5s 
ovemiled and c&Jled At providing a ir&ew shot record. As an exaggerated estampie, fm B 1 
5 secoad sweep dilated by 8 second aiad sampled at a 2 milK&ssoflDd irate (qui? »p|g te&OT 
2 gnilliseeoodls), the sweep is reaampled to 1 mSHiBecosid providing twitee is sa&ay "saaaples. 
Tk iresampkd data is then give?? a sample of 2 omillis^&ot^ds makkg tlhs rasoird ftwks as 
Eeiag. T&s dilated sweep is then correlated She sw>sHiibted sw©gp. The phase of the result 
5s the 5^a3 to3 phase correction. 
XO The ptoa© component of She data is cum-iem due to die Doppfer shi^s which TQBzh 

from the ship's motioa. 'Moreover, the major efifetf of the distortion db© do §m*fi®3 ffltotol as 
ess© ©sJy k the phase cpmrnm of the dsma. Yhz dmmia® can be dammmed by tedng the 
phase ©onBpsas&iat of fthe daaa to a comgaamt vsfee, preferably zero pfoas©. Tte* 5a ascordsaiee 
with She pirefeinred emrabodoment, flte ©state aseaswemem aod processing &ys&©Ea Si @Ktirs6l8 
SS the pfes ©DKrapostesiii of the dilation model cross^cosrelatod with the mefbrerace sweep.. A 
standard all pass Sawsrse filter, suclb as a Wie^^^vinson fiBteff 5s ssJaeGGd to disa5ja®a@ fee 
phase ©osteat of the recorded dste ate eonektoosn with tth© mgferaice sweep signal. Aa all 
pass afflvgsse Ste* doss no* alter tics amplitude coated of the data, r^fe^r ©a8y She pfess 
dttDtsoiL The iter preferably ns cmsgracted to mmv® the phase content, cs^sso^ OTSpw 
20 iigaol tar* <fe SCter 8© tave zsro ptee. 

Ifts <sot©^si ffikm preferably are appifed to the eMire trace of dalfe ©ad (ten tite 
appropriate seginrcsotss from ©acta connected tracts aire sdested m& eoiBsbimed tog^fe* go fossa o 
©s^J^^ gffim^^^!^^ TTIh^ss 9 ftlfcg (^inns^t5®jro for tifes dasa Ifer s eeisrso gas© o>f 1 s^oissdL ife^ 
©ample, h oppfisiD 2© She dosi. SMterOy, the eanrectiGKns tfbir toe geSes «f 2 stands*,. 3 
23 sse^sds, 4 see®^ aad otot* sr© &lso applied to the dteta set, theireby gsiraeimtifljg firar data s^s 
eoefe eors9^©d) by a pixtieylsr correction fslaer. Tte^ only the corre^sd dam 0 oo 1.5 
£sc®ftds h select i5r®m the first tea s@tp the corrected data frooft 8 .5 to 2.S s&wrcds is gelecC@S 
torn secoiod <Sofe set, ths oosnected d@t& from 2.5 to 3 J $$eonds is s^ie^ed from the ahosd 
data s^, ahe ccw^led dfaJa &om 3.5 ao 4.5 seconds is select^ ftonn the fourth dlasa sa, aad so 
30 ©a. 

Afi<3- Goramisig afee seismic dsas for receiver and source motion for a constant dip slice 
m step 210 (IF^Mre 3), the next dip slice of F* data is selected in step 240 sad steps 200-210 
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are repeated until all of the F-K data has been selected. Once all the data has been corrected 
for each dip slice of F-K data, the results are summed in step 230 to produce the desired data 
corrected for source and receiver motion. 

Numerous variations and modifications will become apparent to those skilled in the art 
once the above disclosure is fully appreciated. It is intended that the following claims be 
interpreted to embrace all such variations and modifications. 
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CLAIMS: 

1 . A method for ooiTecting seismic data for distortion caused by the motion of a seismic 
source and the motion of a seismic receiver, comprising: 

(a) correcting the seismic data for motion of the receiver; 
5 (b) correlating the receiver corrected seismic data from step (a) with a reference 

frequency sweep signal used by said seismic source to generate seismic waves; 

(c) computing a forward F-K transform of the correlated data from step (b) to 
produce F-K transformed data; 

(d) selecting a constant dip slice of said F-K transformed data; 

10 (e) computing an inverse F-K transform of said constant dip slice of said F-K 

transformed data to produce inverse F-K transformed data; and 
(f) computing an appropriate correction filter for said inverse F-K transformed data. 

2. The method of claim 1 wherein said step of selecting a constant dip slice of said F-K 
15 transformed data includes selecting F-K transformed data corresponding to a 

predetermined range of dip angles. 

3. The method of claim 2 wherein said steps of selecting a constant dip slice of F-K 
transformed data (step (d)), computing an inverse F-K transform of said selected 

20 constant dip slice of data (step (e)X and computing an appropriate correction filler for 

said inverse F-K transformed data (step (f)) are repeated for additional constant dip 
slices of F-K data. 

4. The method of claim 3 wherein said correction filter is applied to said seismic data to 
25 correct for said source motion. 

5. The method of claim 4 wherein said step of computing an appropriate correction filter 
includes calculating at least one dilation value for said inverse F-K transformed data. 

30 6. The method of claim 5 wherein said seismic data is generated by said source emitting 
seismic waves using said reference frequency sweep signal, and said dilation value is 
calculated by subtracting the time a seismic wave takes to travel from the source to the 
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receiver at the end of the frequency sweep from the time a seismic wave takes to travel 
from the source to the receiver at the beginning of the frequency sweep. 

7. The method of claim 6 wherein said inverse F-K transformed data includes seismic data 
reflected off of a plurality of diffracton located along lines corresponding to said 
predetermined range of dip angles and a plurality of dilation values are calculated for 
said seismic data. 

8. The method of claim 7 wherein said dilation values are discretized into segments. 

9. The method of claim 8 wherein said step of computing an appropriate correction filter 
includes resampling said seismic data to a sampling period that is different than the 
sampling period used to initially sample said seismic data. 

15 10. The method of claim 9 wherein said step of computing an appropriate correction filter 
also includes computing a plurality of all pass inverse filters, each all pass inverse filter 
corresponding to a particular segment of said dilation values. 

11. The method of claim 10 wherein the results of applying said correction filters to said 
20 seismic data are added together. 

12. A method for correcting marine seismic data for distortion caused by the motion of a 
sebniic source, comparing: 

separating said seismic data imo portions, each portion 
25 of dip angles; 

computing Dopplcr correction fitters for said portions; and 
applying said Dopplcr correction filters to said portions 

13. The method of claim 12 wherein said step of separating said seismic data into portions 
30 includes: 

computing a transform of said seismic data to produce transformed data; 
selecting a subset of said transformed data; and 
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computing ar. inverse transform of said subset. 

14. The method of claim 13 wherein said step of computing a transform includes 
computing an F-K transform and said step of computing an inverse transform includes 

5 computing an inverse F-K transform. 

15. The method of claim 14 wherein said step of computing a transform includes 
computing a t-p transform and said step of computing an inverse transform includes 
computing an inverse t-p transform. 

in 

16 The method of claim 15 wherein said step of computing a transform includes 
computing a Laplace transform and said step of computing an inverse transform 
includes computing an inverse Laplace transform. 

IS 17. The method of claim 1 6 wherein said Doppler correction filters include all pass inverse 
filters. 

18. A seismic system, comprising: 

a seismic measurement and processing system; 
20 a user input device coupled to said seismic measurement and processing system; 

a display device coupled to said seismic measurement and processing system; 
wherein said seismic measurement and processing system includes a processing 
unit and a storage unit, and said seismic measurement and processing system corrects 
seismic data for source motion by computing an F-K transform of said seismic data, 
25 computing an inverse F-K transform of a constant dip slice of said F-K transformed 

data, and computing and applying appropriate correction filters for each constant dip 
slice of F-K transformed data. 

19. The seismic system of claim 18 wherein said seismic measurement and processing 
30 system selects a constant dip slice of said F-K transformed data corresponding to a 

predetermined range of dip angles. 
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20. The seismic system of claim 19 wherein said seismic measurement and processing 
system selects multiple pie slices of F-K transformed data, computes a set of 
appropriate correction filters for each selected constant dip slice of data, applies said set 
of appropriate correction filters to said seismic data detected by said receiver, and adds 

5 together the seismic data after being filtered by said set of correction filters. 

21. A seismic system, comprising: 

a seismic measurement and processing system; 

at least one seismic source coupled to said seismic measurement and processing 

ID system; 

at least one seismic receiver coupled to said seismic measurement and 
processing system; 

a user input device coupled to said seismic measurement and processing system; 

and 

13 a display device coupled to said seismic measurement and processing system; 

wherein said seismic measurement and processing system includes a source and 
receiver interface unit, a processing unit, and a storage unit, and said seismic 
measurement and processing system corrects for source motion by computing an F-K 
transform of seismic data detected by said receivers, computing an inverse F-K 

20 transform of a constant dtp slice of said F-K transformed data, and computing and 

applying appropriate correction filters for each pic slice of F-K transformed data. 

22. The seismic system of claim 21 wherein said seismic measurement and processing 
system selects a constant dip slice of said F-K transformed data corresponding to a 

23 predetermined range of dip angles. 

23. The seismic system of claim 22 wherein said seismic measurement and processing 
system selects multiple constant dip slices of F-K transformed data, computes a set of 
appropriate correction filters for each selected constant dip slice of data, applies said set 

30 of appropriate correction filters to said seismic data detected by said receiver, and adds 

together the seismic data after being filtered by said set of correction filters. 
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